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Summary. We present interferometric observations of the Be star 51 Ophiuchi. 
These observations were obtained during the science demonstration phase of the 
MIDI instrument at the Very Large Telescope Interferometer (VLTI). Using MIDI, 
a Michelson 2 beam combiner that operates at the N band (8 to 13 /im), we ob- 
tained for the first time observations of 51 Oph in the mid-infrared at high-angular 
resolution. 

It is currently known that this object presents a circumstellar dust and gas disk 
that shows a very different composition from other Herbig Ae disks. The nature of 
the 51 Oph system is still a mystery to be solved. Does it have a companion? Is it a 
protoplanetary system? We still don't know. Observations with MIDI at the VLTI 
allowed us to reach high-angular resolution (20 mas). We have several uv points that 
allowed us to constrain the disk model. We have modeled 51 Oph visibilities and 
were able to constrain the size and geometry of the 51 Oph circumstellar disk. 



1 Introduction 

51 Ophiuchi is a Be (B9.5IIIe) star, located at 131 pc with a rotational 
velocity vsini — 267 ± 5 km/s [1]. 51 Oph shows a large infrared excess, 
that was first noticed by Waters et al. 1988 [10]. This infrared excess was 
explained as being originated by the circumstellar dust around the star. 

The 51 Oph circumstellar gas and dust origin is still unknown. Since its 
infrared excess was noticed, this object has been observed in many different 
wavelengths, from the infrared to the ultraviolet. In 1993, Grady and Silvis [3] 
found the first evidence for the presence of dust by detecting a 10 /im silicate 
feature in emission. 51 Oph has been compared to (3 Pic due to the presence 
of a gas envelope with a variable column density, high density gas accreting 
to the star and the fact that this gas is also collisionally ionized. There was 
no detection of an extended disk at 18 [im, which suggested that the hot dust 
must be located in the close proximity of the star [4]. Van den Ancker et al. 
[8] analyzed 51 Oph ISO archive data and suggested a few different scenarios 
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for the nature of this system. The authors have studied the composition 
of the circumstellar gas and detected spectral features due to hot gas-phase 
molecules not typical of Ae/Be stars. They have suggested that 51 Oph might 
not be a young system but a highly evolved one instead. Nevertheless, the 
presence of all the circumstellar material can not be explained if the system is 
really an evolved one. The other scenario includes the presence of a companion 
in order to explain the hot gas observed. A more exotic explanation consists 
of assuming that the detected material is a result of a recent collision of two 
gas-rich planets or the accretion of a solid body as the star increases its size 
at the end of its main-sequence life [8] . 

Roberge et al. [7] observed 51 Oph circumstellar disk with the Far Ultra- 
violet Spectroscopic Explorer and found that the composition of the infalling 
gas is highly nonsolar. However, all studies suggest that the circumstellar 
matter is most likely in the form of a Keplerian disk rather than in a spher- 
ical shell. In 2004 the first MIDI observations of Herbig Ae/Be stars were 
published by Leinert et al. [5] . The authors have studied a sample of 7 Her- 
big Ae/Be stars but could not fit a Dullemond et al. (2001) type disk model 
to the 51 Oph data like they did for the other stars. 

In this paper, we have used the data obtained by Leinert et al. [5] on 
51 Oph and constrained the size and geometry of the 51 Oph circumstellar 
disk by using a standard disk model. In section 2 the observations and data 
reduction are presented. The results are introduced and discussed in sections 
3 and 4, and summarized in section 5. 

2 Observations and Data Reduction 

51 Oph was observed during MIDI Science Demonstration Time (SDT) at the 
VLTI on June 15th and 16th 2003. MIDI is a 2 beam combiner that operates 
in the mid-infrared, 8 to 13 /um. The baseline used was the UT1-UT3 (102m), 
with a maximum full spatial resolution of 20 mas at 10 /im. The observations 
log is presented in Table 2. We have obtained 5 visibility points, that were 
dispersed at a resolution of 30 in the N band. 

We have reduced the 51 Oph data using the IDL software written for 
MIDI. The data reduction is described in detail in Leinert et al. (2004) [5]. 
The data analysis and modeling was done with the interpreted programming 
language yorick. 

3 Results 

We have modeled the circumstellar disk by using a flat disk model with a 
temperature law proportional to r~ 3 ' 4 [6]. The disk radial structure was de- 
rived from Lynden-Bcll & Pringle (1974). In this model, the disk is considered 
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Table 1. Log of the observations [?] 



UT Date 


Projected 


P.A. 


Hour Angle 


Calibrators 




Baseline (m] 


(°) 






15/Jun/03 


101.2 


23 


3 


HD 168454, HD 168454, HD 167618 


15/Jun/03 


101.4 


38 


7 


HD 168454, HD 168454, HD 167618 


15/Jun/03 


85.6 


45 


8 


HD 168454, HD 168454, HD 167618 


16/Jun/03 


98.8 


-7 





HD 165135, HD 152786, HD 165135 


16/Jun/03 


99.6 


14 


2 


HD 165135, HD 152786, HD 165135 



Table 2. Disk model parameters for 51 Oph.[?] 



Parameter 


Best fit Value 


Distance 


131 pc 


R* 


7 Rq 


M* 


3.8 M Q 


T c g 


10000 K 


Av 


0.15 


Accretion rate 


7.10" 5 M /yr 


Disk outer radius 


7 AU 


Disk inner radius 


0.55 AU 


Inclination 


88° 


Position Angle 


78° 



as infinitely thin and the flux is calculated through the sum of all the rings 
of radius r and width dr. 

In Table 3 one can find a detailed list of the parameters used to fit the 
standard disk model to our data. The values for the stellar mass (M+), effec- 
tive temperature (T Q s) and atmospheric extinction (Av) were all taken from 
the literature [8]. We have found a good fit for the data with a disk nearly 
edge-on, inclination of 88°, position angle of 78° and an accretion rate of 
7.1CT 5 M /yr. 

For a radius of 7 Rq, we have found a good fit for the spectral energy 
distribution (SED), with an inner disk radius of 0.55 AU and an outer disk 
radius of 7 AU (see Fig. 3). In this figure one can see the SED for 51 Oph 
well reproduced by a standard disk model. The photometry points used in 
these figure were taken from Waters et al. (1988) [10]. 

In Figure 3 on the left, we can see the visibility model for the 51 Oph 
circumstellar disk and the wv-track for our measurements. The 5 visibility 
points are marked with crosses over the wv-track. In this figure north is up 
and west is left, so we note that the position angle is obtained after a Fourier 
transform of the visibility model (-12°+90°=78°). We have also studied the 
visibilities behavior against projected baselines and hour angle. In Figure 3 
on the right, we have plotted the visibilities as a function of hour angle. For 
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Spectral energy distribution 




Wavelength (urn) 

Fig. 1. Spectral energy distribution for 51 Oph. The circles correspond to photom- 
etry measurements compiled from the literature (Waters et al. 1988). The dotted 
line is the stellar contribution, the dashed line is the disk contribution and the solid 
line corresponds to the total energy distribution (star+disk) . Our standard disk 
model successfully reproduces the 51 Oph SED.[?] 




-100 -50 50 100 Hour angle 



Fig. 2. Left:51 Oph disk visibility model. The five visibility points (marked as 
crosses on the uu-track) were obtained in 2 different nights. Right: Visibilities as 
a function of hour angle. Best fit for a standard disk with an inclination of 88° and 
a position angle of 78°.[?] 

both cases we have found a good fit for all the visibility points by using the 
disk model parameters in Table 3. 
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4 Discussion 

The MIDI observations are compatible with a flat circumstellar disk, inclined 
of 88° with a position angle of 78°. These results are compatible with the 
recent results obtained by Thi et al. [9], who have found a high column 
density of CO rotating in a disk in the inner astronomical unit of 51 Oph, 
seen nearly edge-on. We have determined a stellar radius of 7 Rq, a disk 
inner radius of 0.55 AU, an outer radius of 7 AU, a disk inclination angle 
of about 88° and a disk position angle of about 78°. Due to the nature of 
the object, our interpretation relies critically on one visibility point and thus 
drew some criticism. Four of the visibility points are nearly identical, while 
the fifth one is significantly higher. This higher visibility point was obtained 
with the lowest projected baseline at 85.6m. There is consequently a potential 
degeneracy of information between a base length effect and an angle effect. 
In order to secure our data set, we plan to obtain a few extra visibility points 
with MIDI. 



5 Conclusions 

51 Oph was observed for the first time in the mid- infrared at high- angular 
resolution. 

We have modeled 51 Oph visibilities and were able to constrain the size 
and geometry of the 51 Oph circumstellar disk. 

The best fit to the data corresponds to a flat circumstellar disk, inclined 
of 88° with a position angle of 78°. 

We will need more visibility measurements in order to secure our data set 
and confirm the results obtained. 
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